Bidimensional intermittent search processes: an alternative to Levy flights strategies 
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Levy flights are known to be optimal search strategies in the particular case of revisitable targets. 
In the relevant situation of non revisitable targets, we propose an alternative model of bidimensional 
search processes, which explicitly relies on the widely observed intermittent behavior of foraging an- 
imals. We show analytically that intermittent strategies can minimize the search time, and therefore 
do constitute real optimal strategies. We study two representative modes of target detection, and 
determine which features of the search time are robust and do not depend on the specific character- 
istics of detection mechanisms. In particular, both modes lead to a global minimum of the search 
time as a function of the typical times spent in each state, for the same optimal duration of the 
ballistic phase. This last quantity could be a universal feature of bidimensional intermittent search 
strategies. 

PACS numbers: 87.23. -n 



Search processes, involving a searcher and a target of 
unknown position, play an important role in many phys- 
ical, chemical or biological problems. This is for instance 
the case of reactants diffusing in a solvent until they 
get close enough to react or of a protein searching 
for its specific target site on DNA 0j S • One can 
also mention animals searching for food 
or coast-guards trying to locate wreck victims In 
all these examples, it is of great importance to mini- 
mize the search time. Since the pioneering works of 
Viswanathan et al 7], the question of determining op- 
timal search str ateg ies has appealed a growing attention 
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In this context, Levy flights strategies have been 
proved to play a crucial role in such optimization 
problems nil 0. However, two limitations of these 
strategies have to be mentioned. First, Levy flights tra- 
jectories have been shown to optimize the search effi- 
ciency, but only in the particular case where the targets 
are regenerated at the same location after a finite time 
0,01 which can not be taken as a general rule. Indeed, 
in the case of destructive search where each target can 
be found only once, or in the case of a single target, the 
optimal strategy proposed in 7] is not anymore of Levy 
type, but reduces to a linear ballistic motion. Second, as 
for the applications to behavioral ecology, the destructive 
search is relevant to many situations 5', (t|. However, the 
purely ballistic strategy predicted by 7] in that case can 
not account for the generally observed reoriented animal 
trajectories^. 

Alternatively to these Levy strategies, it has been ob- 
served that intermittent search strategies are widely used 
by foraging animals 0, Many searchers combine 

phases of fast displacement, non receptive to the targets, 
and slow reactive search phases. Everyday-life examples 
also confirm that we instinctively adopt such intermittent 
behavior when looking for a lost object: we search care- 



fully around one location, then move quickly to another 
unvisited area and then search again. 

Up to now, only ID models of such intermittent search 
have been developed, providing a satisfactory agreement 
with experimental data from behavioral ecology |inl |. 
Here we develop a model of 2D intermittent search strate- 
gies, which encompasses a much broader field of applica- 
tions, in particular for animal or human searchers. We 
show that bidimensional intermittent search strategies do 
optimize the search time for non revisitable targets. We 
explicitly determine how to share the time between the 
phases of non reactive displacement and of reactive search 
to find a target in the quickest way. From a technical 
point of view, we also obtain as a by-product the mean 
first passage time for a Pearson type random walk, which 
belongs to a class of non trivial problems which have been 
investigated for long El HHHIiil . Our approach 
relies on an approximate analytical solution based on a 
decoupling hypothesis, which proves to reproduce quan- 
titatively our numerical simulations over a wide range of 
parameters. 

Following we consider a two state searcher (see 
Fig^ of position r that performs slow reactive phases 
(denoted 1), randomly interrupted by fast relocating bal- 
listic flights of constant velocity v and random direction 
(phases 2). We assume the duration of each phase i to 
be exponentially distributed with mean r.;. As fast mo- 
tion usually strongly degrades perception abilities |Rll9j. 
we consider that the searcher is able to find a target 
only during reactive phases 1. The detection phase in- 
volves complex biological processes that we do not aim 
at modeling accurately here. However, essentially two 
modes of detection can be put forward, and lead to dis- 
tinguish between two types of reactive phases 1. The first 
one, referred to in the following as the "dynamic mode" , 
corresponds to a diffusive modeling (with diffusion co- 
efficient D) of the search phase as recently proposed in 
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FIG. 1: Two models of intermittent search: The searcher 
ahernates slow reactive phases (regime 1) of mean duration 
Ti, and fast non reactive ballistic phases (regime 2) of mean 
duration T2. Left: the slow reactive phase is diffusive and 
detection is infinitely efficient. Right: the slow reactive phase 
is static and detection takes place with finite rate k. 

[lo| in agreement with observations for vision|23, tac- 
tile sense or olfaction The detection is assumed to 
be infinitely efficient in this mode : a target is found 
as soon as the searcher-target distance is smaller than 
the reaction radius a. On the contrary, in the second 
mode, denoted as the " static mode" , the reaction takes 
place with a finite rate fc, but the searcher is immobile 
during search phases. Note that this description is com- 
monly adopted in reaction-diffusion systemsQ| or oper- 
ational research^^. A more realistic description is ob- 
tained by combining both modes and considering a diffu- 
sive searcher with diffusion coefficient D and finite reac- 
tion rate k. In order to reduce the number of parameters 
and to extract the main features of each mode, we study 
them separately by taking successively the limits k oo 
and D — !■ of this general case. More precisely, in these 
two limiting cases, we address the following questions : 
what is the mean time it takes the searcher to find a tar- 
get? Can this search time be minimized ? And if so for 
which values of the average durations of each phase? 

We now present the basic equations combining the two 
search modes introduced above in the case of a point-like 
target centered in a spherical domain of radius b with 
reflecting boundary. Note that this geometry mimics 
both relevant situations of a single target and of infinitely 
many regularly spaced non revisitable targets. For this 
process, the mean first passage time (MFPT) at a target 
satisfies the following backward equation |22j: 

I r^"" 

D^lh + 7^ / (i2 - h)d9^ - khir)h = -1 (1) 



V-Vri2 {t2-tl)^-l (2) 

where ti stands for the MFPT starting from state 1 at 
position r, and t2 for the MFPT starting from state 2 



at position r with velocity v. Ia(r) = 1 if |r| < a and 
Ia(r) = if |r| > a. In the present form, these integro- 
differential equations do not seem to allow for an exact 
resolution with standard methods. We propose here an 
approximate resolution based on the introduction of the 
following auxiliary functions : 

s(r) = — / tsd^v, d(r) = — / t2^d0^. (3) 

Averaging Eq.Q and Eq.Q times v over 6^, one suc- 
cessively gets 

V.d-l(s(r)-ti) = -Id=^ / (vVtsW^v, (4) 
which gives in turn 

v-d = £Ea5^(-v.).. (5) 

where {■)g^ stands for the average over 6^. We now make 
the following decoupling assumption 

{viVjt2)e^ — {viVj)e^{t2)0^ = Y%s(r) (6) 

which leads, together with Eq.|0J, to the diffusion-like 
equation: 

D\7^.s{r)--{sir)-h)^-l (7) 

T2 

where D = v'^T2/2. Rewriting Eq.QJ as 

DV\ + -(s(r) - h) - kla{r)h = -1, (8) 

Tl 

Eqs.Q and ^ together with vanishing normal deriva- 
tives at |r| — b provide a closed system for the variables 
s and tl, whose resolution is lengthy but standard. The 
validity domain of assumption 10 is much broader than 
the "Brownian" limit v oo and T2 — > with D fixed, 
in which t2 is independent of the direction of v. Indeed, 
it is also valid in the limit vt2 ^ b, in which a ballistic 
phase includes many reorientations due to successive re- 
flections on the boundary r = b. In addition, it can be 
shown that in one dimension this assumption is exact. 

We first present the solution of Eas H7l8|l in the "dy- 
namic mode" (fc — > oo). The search time (i), defined 
as tl uniformly averaged over the initial position of the 
searcher (note that this last averaging reflects the com- 
plete ignorance of the target position), reads in this case: 
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Af L_ a'^Dri (3-4 ln(fe/a))6'* - 40^62 
2L+ L+ sDto - 



(9) 



with i± = la{a/\[DT2) {li(ba)Ki{aa) - Ii(aa)Ki(6a)) ± a\[DT2 Ii(a/\/z)T^) (Ii(6a)Ko(aa) + Io(aa)Ki(&a)) 



n:— aV^^2 r~ 

and M = Io(a/V£'T2) (Ii(6a)Ko(aa) + Io(aa)Ki(6a)) - A—^ —\{a/dDT2) (Ii(6a)Ki(aa) - Ii(aa)Ki(6a)) 
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FIG. 2: Simulations (points) versus analytical approximate 
(line) of the search time in the "dynamic mode": the search 
time rescaled by the value in absence of intermittence as a 
function of T2 (left) and Inri (right) (the logarithmic scale 
has been used due to the flatness of the minimum), for D = 1, 
V = 1, h = 451. Left: a — 10, ri — 1.88 (green boxes); a = 1, 
ri — 50 (blue circles); a — 0.1, ri — 250 (red crosses). Right: 
a — 10, T2 = 21.6 (green boxes); a = 1, T2 — 22.5 (blue 
circles); a = 0.1, T2 = 41.6 (red crosses). 



where a = {1/{Dti) + 1/(Dt2))^^'^, and 1^ and are 
modified Bessel functions. This expression Q has proved 
to be in very good agreement with numerical simulations 
for a wide range of the parameters (see Fig.||2Jl). The 
optimization of the explicit expression Q leads to sim- 
ple forms in the following situations, depending on the 
relative magnitude of the three characteristic lengths of 



the problem a, b, D/v. We limit ourselves to the the case 
of low target density (a ^ 6) , which is the most relevant 
for hidden targets search problems. Three regimes then 
arise. In the first regime a b ^ D/v, the relocating 
phases are not efficient and intermittence is useless. In 
the second regime a D /v <ti b, it can be shown that 
the intermittence can significantly speed up the search 
(typically by a factor 2) , but that it does not change the 
order of magnitude of the search time. On the contrary, 
in the last regime I?/ti <C a <C 6, the optimal strategy, 
obtained for 



D In^ib/a) 
2^21n(6/a) - 1 



r2,„.in^-(ln(V«)-l/2)'/', 

V 

(10) 

leads to a qualitative change of the search time which 
can be rendered arbitrarily smaller than the non inter- 
mittent search time when u — s- cx). This optimal strategy 
corresponds to a scaling law 



Tl,r, 



D 



2, mill 



a2 (2- 



l/ln(&/«))' 



(11) 



which here does not depend on v. 

We now turn to the "static mode" {D — > 0), which 
leads to the following expression for the search time 



where x = 



2fcri a 

1 -I- fcri VT2 



and y = 



2fcri 



1 -I- fcri VT2 



(13) 



Here again, this expression (|12|l is in very good agree- 
ment with numerical simulations for a wide range of the 
parameters (see Fig.©). In that case, intermittence is 
trivially necessary to find the target, and the optimiza- 



tion of the search time H12() leads for b ^ a to 



n, 



a \i/2 /21n(6/a) - 1 



f-V 

\vkJ 



1/4 



T2- 



-(Mb/a)-l/2f\ 



(14) 
(15) 
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FIG. 3: Simulations (points) versus analytical approximate 
(line) of the search time in the "static mode" : the search time 
rescaled by the optimal value as a function of T2/a (left) and 
ri (right), for fc = 1, u = 1, fo = 28. Left: a = 10, ti = 2.14 
(green boxes); a — 1, ti — 0.9 (blue circles); a = 0.1, ri = 0.3 
(red crosses). Right: a — 10, T2 — S (green boxes); a — 1, 
T2 = 1.7 (blue circles); a = 0.1, T2 — 0.22 (red crosses). 

which corresponds to the scaling law T2^inin = 
which stiU does not depend on v. 

The main resuhs (fTn|) and itTH) . ((T^ obtained in the 
two modes of search lead to extract the following remark- 
able characteristics of intermittent search processes: (i) 
In both cases the search time (t) presents a global mini- 
mum for finite values of the r^, which means that inter- 
mittence is an optimal strategy, (ii) A very striking and 
non intuitive feature is that both modes of search lead to 
the same optimal value of r2^min- As this optimal time 
does not depend on the specific characteristics D and k 
of the search mode, it seems to constitute a general prop- 
erty of intermittent search strategies, (iii) The optimal 
Ti.min are different and depend explicitly on D and fc, 
leading to different scaling laws which are susceptible to 
discriminate between the two search modes . 

Finally we remark that this model provides as a by- 
product an approximation for the MFPT for a Pearson 
type random walk in the spherical geometry previously 
defined: the searcher performs ballistic flights reoriented 
at exponentially distributed times, and, as opposed to 
standard Pearson walks, the target can be found only 
when the distance between the target and a reorienta- 
tion point is less than a. This quantity, obtained here 
straightforwardly by taking k —^ oo and ti — > in 
Eg. lfT^ . writes: 

+ -^{b^ Hb/a) + (62 - a2)(a2 _ 352 + 4wV|)^ (16) 

To our knowledge, a similar result for standard Pearson 
walks is still missing. Note that in the limit v 00, T2 
with D ~ v^T2/2 fixed, the approximate expression l(TB|) 
gives back the well known exact expression for the MFPT 
of a Brownian particle between concentric spheres [2^. 
Moreover, for b ^ a, the search time (jlBII is minimized 
again for the same value Hll)|) and H15|) of T2 , in agreement 
with the limit A: ^ 00 of Eqs. ((T3|l . 



To conclude, we have proposed a two state model 
of search processes for non revisitable targets, which 
closely relies on the experimentally observed intermit- 
tent strategies adopted by foraging animals. Using a de- 
coupling approximation numerically validated, we have 
shown analytically that in the physically most relevant 
bidimensional geometry, intermittent strategies minimize 
the search time, and therefore constitute optimal strate- 
gies, as opposed to Levy flights which are optimal only 
for revisitable targets. We studied two representative 
modes of search, and determined which features of the 
corresponding optimal strategies are robust and do not 
depend on the specific characteristics of the search mode. 
In particular both modes lead to a global minimum of 
the search time as a function of the typical times spent 
in each state, and the optimal duration of the ballistic re- 
location phase is the same for these both modes. As this 
last time does not depend on the nature of the search 
mode, it could be a universal feature of bidimensional 
intermittent search strategies. 
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